CPHT-PC 028.0505 



in 
o 
o 

(N 
>v 



Gluon and Ghost Propagators in Landau Gauge on the Lattice 



A.Y. LokhoV 2 , C. Roiesnel 
Centre de Physique Theoriqu^ de I'Ecole Polytechnique, 91128 Palaiseau Cedex, France. 



> 

O 
O 

in 
o 
in 
o 



X 



We study the ultraviolet behaviour of the ghost and gluon propagators in quenched QCD 
using lattice simulations. Extrapolation of the lattice data towards the continuum allows to 
use perturbation theory to extract Aqcd - the fundamental parameter of the pure gauge theory. 
The values obtained from the ghost and gluon propagators are coherent. The result for pure 
gauge SU(3) at three loops is A^g- « 320MeV. However this value does depend strongly upon 
the order of perturbation theory and upon the renormalisation description of the continuum 
propagators. Moreover, this value has been obtained without taking into account possible 
power corrections to the short-distance behaviour of correlation functions. 



1 Introduction 



Correlation functions involving ghosts appear in Schwinger-Dyson equations in covariant gauge. 
Thus, these correlation functions play an important role in the non-perturbative studies of the 
infrared behaviour of the gluon propagator and running coupling constant, and of the confine- 
ment. Lattice simulations provide a non-perturbative way to obtain propagators and vertices 
in Landau gauge. The ghost propagator has not been studied^ so much as the gluon propa- 
gator in the case of QCD (577(3) gauge group). That is why one should first probe the lattice 
formulation of the ghost propagator by studying its behavior at small distances, where lattice 
results may be checked using perturbation theory. Moreover, the perturbative fit of the ghost 
and gluon propagators gives us the value of Aqcd - the fundamental parameter of the pure 
gauge theory, which is a renormalisation group invariant, but is scheme-dependent. Even if this 
parameter is not a physical quantity, it can serve as a definition of the unique energy scale of 
quenched QCD. In this contribution we report on our on-going study of the ghost and gluon 
propagators in the domain of energy 2 GeV <-> 6 GeV. We have fitted these correlation functions 
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using three-loop perturbative expansion in the MOM renormalisation scheme^. This allows us 
to obtain two values for Aqcd; and thus to check the self-consistency of the lattice approach. In 
what follows we will describe briefly the lattice formulation in the case of quenched QCD and 
present preliminary results for two-points correlation functions. 



2 Simulation Setup 

On the lattice one works with an Euclidean formulation of QCD in discretised four-dimensional 
space (with lattice spacing a). In practical simulations one considers this theory in a finite 
volume (V = L ), with periodic boundary conditions for the gluon field. Then one can evaluate 
numerically the functional integrals which define Green functions, using Monte-Carlo techniques. 
In our simulations we have used the standard Wilson action to produce gauge field configurations 
(more details may be found in reviews^. Then the gauge has to be fixed. Presently, the Landau 
gauge is the only covariant gauge which can be implemented efficiently on the lattice. The 
Landau gauge is fixed by finding a minimum of the functional 

Fu({g(x)}) = -ReTr £^ g{x)U^x)g^{x + e M ), g{x) € SU(3), (1) 

where the link variable is related to the gauge field A®(x) by U^x) = exp zA^(x)t a , t a are the 
generators of the gauge group. A local minimum defines explicitly the gauge transformation 
that transforms {U} to {U^} so, that uff'(x) satisfies the discretised Landau condition 

J ±^F u ({g(x)})=0 ReTr^2[ujf\x) + U^l(x-e,)] = 0. (2) 

Once the gauge is fixed, one can calculate numerically the gluon two-point correlation function 
< A*(x)A b ^x) >= || [VU]Al(x)A b ^y)e- s ^ u \ (3) 

where A^(x) = ^mO") _ j n momentum space, and in the continuum limit, the general 

parametrisation of the gluon propagator in Landau gauge reads 

< A^p)Al(-p) >= ^V(<W - ^). (4) 

Thus, lattice simulations can provide non-perturbative information about the scalar function 
G( P 2 ). 

The ghost propagator in covariant gauge is defined in the continum formulation of the QCD 
as the inverse of Faddeev-Popov operator 

M ab (x,y) = [d-D] ab (x,y). (5) 

where D is the covariant derivative in the adjoint represnetation. On the lattice this definition 
should be modified. Indeed, the Faddeev-Popov operator on the lattice is defined by the second 
derivative of the functional Fjj at a local minimum 

J^^Fv (W«)» -*S*.C«.V). W 
where (f abc are the structure constants of the gauge group) 

M&teix, y) = J2{G^(x){S y>x+fi - S y>x ) + G ab {x - A)(<W " <W) 

+ \ E f abc {5y, x+ »Afa) - dy^A^x - fx)) | 



(7) 
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Figure 1: Three-loop fit for the data obtained for V 
[2, 5GeV; 4GeV] for the ghion propagator (left image), ar 
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; 32 2 ,/3 = 6.4. The fit has been done in the interval 
[2, 5GeV; 4GeV] for the ghost propagator (right image). 



and 

Gf(x) = l -T? ({f ,t>} (U,{x) + Ul(x))) (8) 
Inverting this operator with a local source at the origin with no zero-mode component 

l-I-I -- 

(V is the volume of the lattice), and performing a Fourier transformation, one can extract the 
scalar form factor F(p 2 ) of the ghost propagator in Landau gauge for all momenta at once 

< c a {p)c\-p) >= ^fV. (9) 

Let us discuss briefly the errors of the method. Obviously, one has the statistical error due to 
the Monte-Carlo calculation of the functional integral. The systematic error is mostly due to the 
space-time discretisation. Indeed, the rotational symmetry is broken down to the group H4 - the 
symmetry group of the hypercube. That means that, at most momenta, the calculated values 
of any scalar quantity will split, because an orbit of the H4 group is not uniquely defined by 
the value of p 2 . Thus, one should extrapolate all data points to the continum limit a — > 0. Now 
we have a well developed method of performing such an extrapolation^^] All plots, presented 
further, have been done using this method. Another source for systematic error comes from the 
Euclidean formulation of the QCD and gauge fixing (lattice version of Gribov copies^ ). But 
the so-called Gribov noise is negligible in the short-distace study that is discussed here. 



3 Results 



We have done simulations on lattices with size {16 4 , 24 4 , 32 4 } (1000,500,250 configurations, 
respectively) and for the respective values of the lattice bare couplings (3 = {6.0, 6.2, 6.4} [3 = 

where go is the value of the bare coupling at the cut-off energy scale a -1 ). Thus, these lattices 
have roughly the same physical volume. We did a fit of both propagators at three loops in the 
MOM scheme^, which was shown^to provide a good description of the lattice gluon propagator, 



and then converted the value of A,-——- to the more familiar MS scheme. The results for Attq 

MOM 

are given in Tabled Two errors are quoted. The error between parentheses is the statistical one 
and is rather small. The other error is the systematic error which can be quite large, especially 
for the ghost propagator where it can reach 10%. The systematic error has been estimated from 
the stability of the fits in the region [2GeV; 6GeV]. Examples of the fits are given in Fig. ^ 



Table 1: Results of the fits for Aqcd 



p 


L 


a (3)gluon 
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X 2 /d.o.f 


A (3)ghost 
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rVd.o.f 


6.0 


16 


324(2)+^ MeV 


1.2 


322(8)i^MeV 


0.5 


6.2 


24 


320(2)if 4 MeV 


0.9 


326(5)+^MeV 


0.6 


6.4 


32 


312(l)+^ 5 MeV 


1.4 


331(4)+^MeV 


0.9 



4 Conclusion 

We obtain coherent results for Aqcd at three loops from the ghost and gluon propagators. 
This means that the lattice QCD approach is self-consistent in the ultraviolet domain when 
studying correlation functions which involve ghosts. The value for Aqcd that we report here 
is approximatively equal to 320MeV in the MS scheme. However it should be emphasized that 
this value does depend strongly upon the order of perturbation theory used in the fits. This 
dependence manifests itself in two ways. First, two-loop perturbation theory is also able to 
reproduce the data but with a much higher Aqcd scale. Secondly, at any given order, the value 
of Aqcd relies very much on the renormalization scheme used to describe the propagators in the 
continuum. The value of the A^jg scale that we quote here should therefore be considered as an 
effective Aqcd scale. 

In order to check the asymptoticity and the coherence of Aqcd more precisely one should 
also study the ultraviolet behaviour of the gluon and ghost propagators at four loops in several 
MOM schemes. Possible power corrections^ could also make the value of the Aqcd scale lower. 
We currently work on this, and the results will be reported elsewhere. 

Acknowledgements 

We thank our colleagues of the LPT-Orsay, P. Boucaud, J. P. Leroy, A. Le Yaouanc, J. Micheli 
and O. Pene, for many useful discussions. 

References 

1. H. Suman and K. Schilling, Phys. Lett. B 373 (1996) 314 |arXiv:hep-lat /9512003| . 

A. Sternbeck, E. M. Ilgenfritz, M. Muller-Preussker and A. Schiller, arXiv:hep-lat/0412011 

2. I. Montvay, G. Munster, "Quantum Fields on a Lattice," Cambridge Univ. Press (1994). 

3. D. Becirevic, P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodriguez-Quintero and 
C. Roiesnel, Phys. Rev. D 60 (1999) 094509 |arXiv:hep-ph /9903364l. 

4. D. Becirevic, P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodriguez-Quintero and 
C. Roiesnel, Phys. Rev. D 61 (2000) 114508 |arXiv:hep-ph/9910204| . 

5. T. D. Bakeev, E. M. Ilgenfritz, V. K. Mitrjushkin and M. Mueller-Preussker, Phys. Rev. 
D 69 (2004) 074507 |arXiv:hep-lat/0311041 | . 

6. P. Boucaud et al, JHEP 0004 (2000) 006 |arXiv:hep-ph/0003020| . 

7. K. G. Chetyrkin and A. Retey, "Three-loop three-linear vertices and four-loop MOM beta 
functions in massless QCD," arXiv:hep-ph/0007088 



